CONTINUOUS-TIME QUANTUM WALKS AND TRAPPING 
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Recent findings suggest that processes such as the electronic energy transfer through the pho- 
tosynthetic antenna display quantal features, aspects known from the dynamics of charge carriers 
along polymer backbones. Hence, in modeling energy transfer one has to leave the classical, master- 
equation- type formalism and advance towards an increasingly quantum-mechanical picture, while 
still retaining a local description of the complex network of molecules involved in the transport, say 
through a tight-binding approach. 

Interestingly, the continuous time random walk (CTRW) picture, widely employed in describ- 
ing transport in random environments, can be mathematically reformulated to yield a quantum- 
mechanical Hamiltonian of tight-binding type; the procedure uses the mathematical analogies be- 
tween time-evolution operators in statistical and in quantum mechanics: The result are continuous- 
time quantum walks (CTQWs). However, beyond these formal analogies, CTRWs and CTQWs 
display vastly different physical properties. In particular, here we focus on trapping processes on a 
ring and show, both analytically and numerically, that distinct configurations of traps (ranging from 
periodical to random) yield strongly different behaviours for the quantal mean survival probability, 
while classically (under ordered conditions) we always find an exponential decay at long times. 
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I. INTRODUCTION 



The transport properties of excitons in organic as well 
as in inorganic molecular solids are of fundamental inter- 
est [Kempe, 2003; Woerner et ai, 2004; Miilken et ai, 
2006; Sillanpaa et ai, 2007; Olaya-Castro et ai, 2008; 
Zhou et ai, 2008]. In general, at high temperatures the 
transport is incoherent and can be efficiently modeled 
by continuous-time random walks (CTRWs) over sets 
of participating centers (atoms, molecules, etc.) [Mon- 
troU & Weiss, 1965]. In this case the transport follows a 
master equation. The transfer rates between the partici- 
pating centers can be related to the spatial arrangement 
of the centers. The arrangement is captured by the so- 
called Laplacian Matrix L which we will identify here 
with the transfer matrix of the CTRW. However, when 
dealing with quantum particles at low densities and low 
temperatures, decoherence can be suppressed to a large 
extent: The study of transport in this regime requires 
different modeling tools, able to mimic the coherent fea- 
tures. Clearly, quantum mechanical transport phenom- 
ena follow Schrodinger's equation. In order to make con- 
tact to CTRW we relate the Hamiltonian of the system 
to the classical transfer matrix, H = — L; this yields a 
description mathematically closely connected to the clas- 
sical master equation approach. Indeed, this realizes a 
quantum mechanical analog of the CTRWs defined on 
a discrete structure, i.e. the so-called continuous-time 
quantum walks (CTQWs). However, apart from formal 
analogies, coherence can give rise to very peculiar prop- 
erties (e.g. Anderson localization [Anderson, 1958], cru- 
cial dependences of the transport on the starting posi- 
tion [Miilken et ai, 2006] and a quadratic speed-up of 
the chemical distance covered [Agliari et ai, 2008]) with 



no counterpart in classical transport. These effects al- 
low interesting and cross-disciplinary applications and 
can also be exploited in experiments in order to distin- 
guish whether the transport is rather classical or rather 
quantum mechanical [Miilken et ai, 2007]. 

In particular, a common means for probing the trans- 
port relies on the interaction with other species, such as 
impurity atoms or molecules (found in or doped into the 
medium) which irreversibly trap the charges or quench 
the excitations. Consequently, a great deal of recent the- 
oretical work has focused on investigating essential fea- 
tures of basic trapping models, wherein a single particle 
moves in a medium containing different arrangements of 
traps. Indeed, much is known about the decay when the 
motion is incoherent [Van Kampen, 1981; Blumcn et al, 
1983; 1986; ben-Avraham & Havhn, 2000], while (as we 
will show here) when quantum effects become important, 
strong deviations from the classical results occur. 

In a set of early works the dynamics of coherent exci- 
tations on a chain with randomly distributed traps has 
been investigated using several methods [Hemenger et ai, 
1974; Kenkre, 1978; Huber, 1980; Parris, 1991] which 
provided a reasonably description of the process at short 
times and in the asymptotic limit. On the other hand, 
from the experimental point of view, the most relevant 
regime is the one of intermediate times; in this time inter- 
val some striking effects have been recently highlighted 
[Miilken et ai, 2007; Miilken et ai, 2008]. 

Here we focus on trapping processes taking place on a 
finite ring where the traps are distributed according to 
different arrangements: the traps are either gathered in a 
cluster or distributed periodically or randomly. In these 
cases the classical survival probability Pmit) has been 
studied intensively (see e.g. [ben-Avraham & Havlin 
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2000]). In fact, under ordered conditions PM{t) is known 
to exponentially decay to zero. Conversely, for a ran- 
dom distribution of traps Puit) exhibits different time 
regimes; At long times it follows a stretched exponen- 
tial which turns into a pure exponential when finite size 
effects dominate. As for the CTQW, the emergence of 
intrinsic quantum-mechanical effects, such as tunneling, 
prevents the decomposition of the problem into a col- 
lection of disconnected intervals and, as we will see, the 
mean survival probability TiM{t) is strongly affected by 
the trap arrangement. Hence, by following the tempo- 
ral decay of IIm {t) we can extract information about the 
geometry. Furthermore, we show that in the cases an- 
alyzed here Pmit) and IlMit) exhibit qualitatively dif- 
ferent behaviours; this allows to determine the nature, 
either rather coherent or rather incoherent, of the trans- 
port process. 

Our paper is structured as follows: In Sec. [TTlwe pro- 
vide a brief summary of the main concepts and of the 
formulae concerning CTQWs. In Sec. IIIII we introduce a 
mathematical formalism useful for analyzing trapping in 
the CTQW picture. In the following Sec. lIV| we consider 
special arrangements of traps on a ring and we investigate 
the mean survival probability by means of a perturbative 
approach; these analytical findings are corroborated by 
numerical results. In Sec.|V]we study the case of random 
distributions of traps and finally, in Sec. I VII we present 
our comments and conclusions. 



II. CONTINUOUS TIME QUANTUM WALK 

Let us consider a graph Q made up of N nodes and al- 
gebraically described by the so-called adjacency matrix 
A:The non-diagonal elements Aij equal 1 if nodes i and 
j are connected by a bond and otherwise; the diago- 
nal elements An are 0. From the adjacency matrix we 
can directly derive some interesting quantities concerning 
the corresponding graph. For instance the coordination 
number of a node i is Zi — Aij and the number of 
walks of length £ from i to j is given by {A^)ij [Biggs, 
1974]. 

We also define the Laplacian operator L according to 



-Aij] the set of all N eigenvalues of L is called 
the Laplacian spectrum of G- Interestingly, the Laplacian 
spectrum is intimately related not only to dynamical pro- 
cesses involving particles moving on the graph, but also 
to dynamical processes involving the network itself; these 
include energy transfer and diffusion-reaction processes 
as well as the relaxation of polymer networks, just to 
name a few (see for example [Mohar, 1991; Galiceanu & 
Blumen, 2007] and references therein). 

In the context of coherent and incoherent transport 
it is worth underlining that, being symmetric and non- 
negative definite, L can generate both a probability con- 
serving Markov process and a unitary process [Childs & 
Goldstone, 2004; Miilken & Blumen, 2005; Voha et at, 
2006]. 



Now, continuous-time random walks (CTRWs) [Mon- 
troU & Weiss, 1965] are described by the following Master 
Equation: 



N 



dt 



(1) 



being Pk,j{t) the conditional probability that the walker 
is on node k when it started from node j at time 0. If 
the walk is unbiased the transmission rates 7 are bond- 
independent and the transfer matrix T is related to the 
Laplacian operator through T = — 7L (in the following 
we set 7=1). 

We now define the quantum-mechanical analog of the 
CTRW, i.e. the CTQW, by identifying the Hamihonian 
of the system with the classical transfer matrix, H = 
-T [Farhi & Gutmann, 1998; Miilken & Blumen, 2005]. 
Hence, given the orthonormal basis set ]j), representing 
the walker localized at the node j, we can write 



N 



N 



H 



(2) 



j=l fcNNj 



where in the second sum k runs over all nearest neighbors 
(NN) of j. The operator in Eq. [2]is also known as tight- 
binding Hamiltonian. Actually, the choice of the Hamil- 
tonian H is, in general, not unique [Childs & Goldstone, 
2004] and Eq. [2]has two important advantages: It allows 
to take into account the local properties of the arbitrary 
substrate and, remarkably, it yields a mathematical for- 
mulation displaying important analogies with the classi- 
cal picture. In fact, the dynamics of the CTQW can be 
described by the transition amplitude akj (t) from state 
\j) to state ]fc), which obeys the following Schrodinger 
equation: 



N 

1=1 



(3) 



structurally very similar to Eq. [TJ The solution of Eq. [3] 
can be formally written as 



ak,jit) = (fclexp(-^m)lj), 



(4) 



whose squared magnitude provides the quantum mechan- 
ical transition probability iTkjit) = \ak,j{t)\'^. 

In general, it is convenient to introduce the orthonor- 
mal basis ]$„), n € [1, N] which diagonalizes H; the cor- 
respondent set of eigenvalues is denoted by {£'n}ri=i,...,Ar. 
Thus, we can write 



N 

E( 

n=l 



|e-^"*l<i>„)($„lj) 



(5) 



It should be underlined that while both problems 
(CTRW and CTQW) are linear, and thus many results 
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obtained in solving CTRWs (eigenvalues and eigenfunc- 
tions) can be readily reutilized for CTQWs, the phys- 
ically relevant properties of the two cases differ vastly: 
Thus, in the absence of traps CTQWs are symmetric 
under time-inversion, which precludes them from at- 
taining equipartition for the TTkjit) (such as the Pk,j(t) 
for CTRWs) at long times. Also, the quantal system 
keeps memory of the initial conditions, exemplified by 
the occurrence of quasi-revivals [Miilken & Blumen, 2005; 
Miilken & Blumen, 2006]. 



III. CTQWS IN THE PRESENCE OF TRAPS 

As discussed in the previous section, the operators de- 
scribing the dynamics of CTQWs and of CTRWs share 
the same set of eigenvalues and of eigenstates. However, 
when new contributions (arising e.g. from the interac- 
tion with external fields or absorbing sites) are incorpo- 
rated, the eigenvalues and the eigestates start to differ. 
In the following we introduce a formalism useful to ana- 
lyze the dynamics of CTQWs and CTRWs in the pres- 
ence of traps; for this we will denote with Hq and Tq the 
unperturbed operators without traps. 

Let us consider a system where M out of the N nodes 
are traps; we label the trap positions with rrij, with j — 
1, ...,M, and we denote this set with Ai. 

For substitutional traps the system can be described by 
the following effective (but non-Hermitian) Hamiltonian 
[Miilken et al, 2007] 



H = Ho iV, 
where T is the trapping operator defined as 



M 



(6) 



(7) 



The capture strength Tm, determines the rate of decay 
for a particle located at trap site m-,-; here we will take 
the Tm- to be equal for all traps, i.e. Tm = T for all 
j. The limit F ^ oo corresponds to perfect traps, which 
means that a classical particle is immediately absorbed 
when reaching any trap. 

Due to the non-hermiticity of H, its eigenvalues are 
complex and can be written as — ei — iji {I — 1, N); 
moreover, the set of its left and right eigenvectors, |$/) 
and ($; I , respectively, can be chosen to be biorthonormal 
((<i>;|$J) = dij') and to satisfy the completeness relation 
E^i = 1- Therefore, according to Eq. O the 

transition amplitude can be evaluated as 

N 

1=1 

from which TTk,j{t) = |Q;/c.j(i)p follows. 



Of particular interest, due to its relation to experimen- 
tal observables, is the mean survival probability Hixiit) 
which can be expressed as [Miilken et ai, 2007] 



n 



M 




N-M ^ 

1=1 



N- 



N 



M ^ 

Ll' = l 



|),|TO)(m|$;) 



{^i,\m){m\^i 



2 

1(9) 



The temporal decay of IIm (t) is determined by the imag- 
inary parts of El, i.e. by the ji. As shown in [Miilken et 
ai, 2007] at intermediate and long times and for M <^ N 
the T\M{t) can be approximated by a sum of exponen- 
tially decaying terms: 



n 



1 



N-M 



N 

1=1 



(10) 



and is dominated asymptotically by the smallest 7/ val- 
ues. 

Now, in the incoherent, classical transport case trap- 
ping is incorporated into the CTRW according to 



= J-o 



(11) 



The transfer operator T is therefore real and symmetric, 
and it leads to real, strictly negative eigenvalues which 
we denote by —A;; to them correspond the eigenstates 

Analogously, the mean survival probability for the 
CTRW can be written as 



PM{t) 



E E p'^j W 



1 



N 



M ^ 

1=1 



N 



(12) 



From Eq. [T2]one may deduce that PM{t) attains in gen- 
eral rather quickly an exponential form; furthermore, if 
the smallest eigenvalue Amin is well separated from the 
next closest eigenvalue, Pmit) is dominated by Amin and 
by the corresponding eigenstate |<^min) [Miilken et al., 
2007; Miilken et ai, 2008]: 



PMit) 



1 



N ~M 



E (^i'^" 



(13) 



Lower estimates of the gap A between the two smallest 
eigenvalues have been found in the past for special choices 
of operators (see e.g. [Chen M., 1997] and references 
therein). For instance, the operator Tq has Amin = 0; 
its next smallest eigenvalue represents the algebraic con- 
nectivity of the graph, namely the relative number of 
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edges needed to be deleted to generate a bipartition. In 
the case of a fc-regular graph A is bounded from be- 
low by k/{DN), being D the diameter of the graph, i.e. 
the maximum distance between any two vertices [Chung 
F.R.K., 1996]. 



independently of the trap arrangement and 

Vi,N-i = -ir/Nj:^i, exp{2^7^mJ[^ - {N - l)]/N} 

= -iV/NY.f^^ey.Y'i'^iTTlmj/N). (21) 
By inserting the last results into Eq. [12] we get 



IV. PERTURBATIVE APPROACH FOR 
TRAPPING ON A RING 

When the strength F of the trap is small with respect 
to the couplings between neighbouring nodes (which here 
means T <C 1), we can treat the effective Hamiltonian 
introduced in Eq. [6] along the lines of time-independent 
perturbation theory. 

Before developing this strategy we fix the structure 
by considering a ring of length N so that the coordination 
number equals 2 for all sites (Z = 21), where we assume 
N to be even. For the corresponding Hamiltonian Ho we 
know exactly all the eigenvalues and eigenvectors; one has 
namely 



e[^^ =2-2cos(2ttI/N) 



and 



(0)n 



1 ^ 



-i2TTlj/N 



(14) 



(15) 



We underline that all the eigenvalues, apart from Ej^ ^2 — 
4 and En — 0, are two-fold degenerate, Ei — E^-i {I — 
1,2,..., N/2 — 1). We now apply perturbation theory to 

evaluate to first order the corrections e[^^ to the eigen- 



values El. For / = N/2 and for / 
degenerate expression 



I 

N we use the non- 



E, 



(1) 



(rn|$ 



(0)n 



(16) 



M± 



M 



(22) 



(1) 



We notice that for special trap arrangements the 
can be calculated exactly: The most striking re- 
sults are obtained when the exponential in the sum 
in Eq. [21] equals one of the values from the set 
{l,i,— 1,— i}. Then the absolute value of the sum re- 
duces to I ^^-^j^ exp(i47rZmj/iV)| = M. For this there 
have to exist indices I 7^ N/2 and / 7^ N such that rrij 
can be expressed as 



N 



(23) 



where kj and c are arbitrary integers and r = 0, 1, 2 or 
3, corresponding to or —i, respectively. Conse- 

quently, we obtain for the correction 



= -^r| and E^^l, = 0, 



(24) 



so that the degeneracy is always lifted. 




and get 



E 



N/2 



M 

4 - iF — and En 

N 



iV- 



N' 



For / different from N/2 and N we set 



(17) 




(0)\ 



(18) 



and we apply the expression valid for two-fold degenerate 
solutions of Ho : 



eI^'^ - 2 + ^N-LN-l) (19) 



1 



±- 



(Vu-VN-LN-lf + m.N-l? 



1/2 



where we choose the positive sign for I £ [1, N/2 — 1] and 
the negative sign for I E [N/2 + 1, — 1]. Now we have 



Vii = Vn-ln-i 



-iV- 



TV' 



(20) 



FIG. 1: Examples of periodic (top) and sequential (bottom) 
arrangements of M = 5 traps on a ring of size A'^ = 20. 

Let us now focus on a periodic distribution of traps 
with mj = jN/M, while N/M S N. It is easy to see that 
in this case there exists a non-empty set T of distinct 
values of I G [1, A^/2— 1] satisfying the condition of Ea.[23l 
this occurs for 21 /M G N, so that the cardinality of T is 
given by the number of integers in {2//M}/^i 2,...,JV/2-i) 
namely by 



ITI = 



L(iV-2)/MJ 
[{N - 2)/2M\ 



for even M, 
for odd M, 



(25) 
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where lx\ denotes the largest integer less than or equal 
to X. In particular, for both M = 1 and M = 2 we have 
|T| = N/2—1. Now, the numerical diagonalization of the 
Hamiltonian H shows that for I G T we get jn-i — (not 
only in first order in T). Consequently, the corresponding 
term in Eq. [9] decays to a non vanishing value, and from 
Eq. [TO] we have for t cx) : 



be written as 



11m (t) 



N- M 



(26) 



Hence, recalling Eq. [551 for large structures with M ^ 
N, IlM{t) asymptotically decays to 1/M (even case) and 
to 1/(2M) (odd case). Figure [51 shows results obtained 
for a ring of size N — 300 with a periodic arrangement 
of M = 10 (|T| 29) and M = 75 (|T| = 1) traps. 
Consequently, the survival probability nj\/(i) decays to 
the constant values 1/10 and 1/225, respectively. From 
a physical point of view, the finite limit for the survival 
probability stems from the existence of stationary states 
to which the nodes in Ai do not contribute, so that they 
never "see" the traps. This genuine quantum-mechanical 
effect has no counterpart in the classical case where, for 
finite structures, the survival probability always decays 
to zero in the presence of traps. In particular, as shown 
in Fig.[2l PM{t) decays exponentially, as expected. 




500 1000 

t 



1500 



FIG. 2: Survival probabilities IlM{t) (continuous line) and 
PM(t) (dotted lines) on a ring of size TV = 300 and F = 0.01 
in the presence of M = 10 and of M — 75 traps arranged 
periodically, i.e. rrij — jN/M. Note the semilogarithmic 
scales. 



Let us now focus on another special configuration of 
M traps, M > 1: we consider a sequential arrangement, 
such that nij — j and j = 1,....,M. Hence, Eg. [5T] can 



M 



Vi^N-i = -ir/N"^exp{Ainlj/N) 
-iF exp(4i7TlM/N) - 1 



(27) 



N exp{4TTil/N) - 1 
-iT s\tl{2tiMI/N) 



exp(47riZ/TV) 



N sin(27r//Af) 
from which we get 



exp[2i7r;(M+ 1)/A^], 



(1) 



_ / sin(27rAf//jV) \ 
TV I sm{2Til/N) ) 



(28) 



We notice that since I ^ N/2 and I ^ N then 21 /N ^ N, 



while for 21M/N e N then e\'^^ = E)^'_i = -iVM/N. In 
particular, when M = N/2, we have 7; = M/N for each 
value of I e [1, A^]- As a result, in Eq. [Sjthe first term 
vanishes due to the completeness property and the fact 
that the 7; are no longer /-dependent. As for the second 
term, by neglecting oscillations, we get 



^ ' N ~M 2 



Tt 



(29) 



which is independent of TV. As shown in Fig. [3| the ex- 
ponential behaviour predicted by Eq. [511 holds also for 
intermediate times. 

In Fig. [31 we compare the survival probabilities of 
CTQWs and CTRWs: as highlighted by the semi- 
logarithmic plot, the decay is exponential in both cases, 
although faster in the former. Indeed, for the CTRWs 
we have in the long-time limit from Eq. 1121 



M \ ) ^ 2 



■rt/2 



(30) 



where we used the fact that the smallest eigenvalue is 
TM/N. By comparing Eq. [29| and Eq. [30] we see that, 
although the decay is exponential in both cases, the decay 
rate is twice larger for nj\/(t) than for PM{t). 



RANDOM DISTRIBUTIONS OF TRAPS 



We now take N to be odd (so as not to fulfill Eq. 
and consider random arrangements of traps: we pick 
M distinct trap locations randomly from a uniform dis- 
tribution and determine the corresponding UMit) and 
P]\i{t). Then we average these over different, indepen- 
dent realizations to determine (Ha/ (i)) and (P/\/(t)). As 
already mentioned in Sec. Ill {PAiit)) exhibits different 
behaviours: in an infinite system the decay law is a 
stretched exponential at long times, whereas in finite sys- 
tems at such times the decay gets to be exponential. In 
Fig. [5] we show evidence of the long-time exponential be- 
haviour of {PM{t)) in systems of relatively small size. 

Let us now consider (J\M{i)) for random trap arrange- 
ments. Now, the (J\M{t)) decay differs qualitatively from 
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1000 2000 3000 4000 

t 



FIG. 3: Survival probability IlM{t) on rings of size A'^ — 
32,48,64 and 96 for T = 0.04,0.01,0.004, as indicated. The 
number of traps is M = N/2 and they are placed consecu- 
tively, i.e. m.j = j. The straight lines represent Eg. 1291 




FIG. 4: Survival probabilities IlM{t) and Phiit) on rings of 
size iV = 48 for r = 0.001; the number of traps is M = N/2 
and they are placed consecutively, i.e. rrij = j. The straight 
lines represent Eq. [5^ (continuous line) and Eq. [30] (dashed 
line). 

that of the Hm (t) analyzed in the previous section. As 
shown in Fig. [HI for intermediate times the average sur- 
vival probability displays a power law, which decays more 
slowly than exponentially: 

(IlMit)) --t-'/^. (31) 

A similar result has already been obtained for CTQWs on 
a finite chain with two traps at its ends {Ai — {l,iV}), 
in the presence of either nearest-neighbour [Miilken et 
ai, 2007] or long-range interactions [Miilken et ai, 2008]. 




FIG. 5: Average survival probabilities {PM{t)) for rings of 
sizes A'' = 51 and 101. Here P = 0.1 and M is either 4 or 
8. The data presented have been averaged over 120 different 
realizations, see text. The straight lines highlight the expo- 
nential decay. 

There one could understand the power law decay based 
on the imaginary part of the Hamiltonian spectrum {7;}, 
which in a large interval scales algebraically with By 
fitting the numerical data obtained for different sizes and 
concentrations we get the characteristic exponent /x de- 
picted in the inset of Fig. [6l 



VI. CONCLUSIONS 

In conclusion, we have modeled the coherent dynam- 
ics by continuous-time quantum walks and studied in- 
teractions with traps: Taking a periodic chain as sub- 
strate, we calculated the mean quantal survival proba- 
bility nA/(i) and we compared it to the classical Pm(0 
for different trap arrangements. The quantum problem 
was approached both analytically (by means of pertur- 
bative theory) and numerically, showing that the spa- 
tial distribution of the traps significantly affects IIm (t) ■ 
In particular, when the traps are arranged periodically 
throughout the substrate, Il]\i(t) decays asymptotically 
to a nonvanishing value which depends directly on the 
system size N and on the number of traps M (e.g., when 
M = 2, n2{t) 1/2 for t 00). This is a genuine 
quantum-mechanical effect with no counterpart in clas- 
sical mechanics, where PM{t) decays to zero for finite 
systems. 

Another interesting, deterministic trap configuration is 
realized by distributing the traps consecutively such to 
form a cluster; then at intermediate and long times the 
survival probability decays exponentially with the char- 
acteristic time r~^. Now, for the same trap configu- 
ration, the characteristic time for the classical survival 
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t 



FIG. 6: (nM(t)) and jj, for rings of sizes A'' = 51, 101, and 201. Here F = 0.1 and M is 2, 4, or 8, as indicated. The data presented 
have been averaged over 120 different realizations. The main figure displays the average survival probabilities (nM(i)) in double 
logarithmic scales. The straight lines represent the best fit. The inset displays the exponent as a function of c = M/N for 
systems of sizes N = 201 (*) and N = 101 (•). 



probability doubles, being 2r^^. 

When the traps are distributed randomly on the sub- 
strate, a further, qualitatively different behaviour of 
IlM{t) is obtained. In fact, by averaging over differ- 
ent independent configurations we find in this case that 
at intermediate times (IlM{i)) decays algebraically, i.e. 
lJlM{t)) ^ t^^/^, where ^ depends on M and N and is 
related to the imaginary part of the Hamiltonian spec- 
trum. On the other hand, for systems of relatively small 
size we find that in the same time range finite-size effects 
dominate {PM(t)), giving rise to an exponential decay. 

These results establish that studying the decay due to 
trapping is indeed an advantageous means to monitor the 
system's evolution, as it allows to determine the nature 
of the transport, which can be either rather coherent or 
rather incoherent. Moreover, the behaviour exhibited by 
njv/(t), being qualitatively affected by the trap configu- 
rations, may be used to distinguish between these. 
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